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The hairpin ribozyme is a prototype small, self-cleaving RNA motif. It exists naturally as a four-way RNA junction containing 
two internal loops on adjoining arms. These two loops interact in a cation-driven docking step prior to chemical catalysis 
to form a tightly integrated structure, with dramatic changes occurring in the conformation of each loop upon docking. 
We investigate the thermodynamics and kinetics of the docking process using constructs in which loop A and loop B 
reside on separate molecules. Using a novel CD difference assay to isolate the effects of metal ions linked to domain 
docking, we find the intermolecular docking process to be driven by sub-millimolar concentrations of the exchange- 
inert Co(NH 3 ) 6 3+ . RNA self-cleavage requires binding of lower-affinity ions with greater apparent cooperativity than 
the docking process itself, implying that, even in the absence of direct coordination to RNA, metal ions play a catalytic 
role in hairpin ribozyme function beyond simply driving loop-loop docking. Surface plasmon resonance assays reveal 
remarkably slow molecular association, given the relatively tight loop-loop interaction. This observation is consistent 
with a "double conformational capture" model in which only collisions between loop A and loop B molecules that are 
simultaneously in minor, docking-competent conformations are productive for binding. 



Introduction 

Due to the polyanionic nature of RNA, the nature and concen- 
trations of cations present in solution play a major role in deter- 
mining the structural and functional properties of a given RNA 
system. 1 " 3 Biophysical analyses of cation binding to specific sites 
in RNA, however, are greatly complicated by the much larger 
numbers of cations interacting in non-specific or atmospheric 
fashion with the phosphodiester backbone. 4 Metal-dependent 
catalytic RNA molecules, or ribozymes, can provide an impor- 
tant window into cation-RNA interactions since the catalytic 
activity can be used as a readout for cation binding and the for- 
mation of appropriate RNA structure. The investigation of the 
variety of structural and functional roles played by metal ions 
in ribozyme mechanisms has been a fruitful source of insight 
into RNA structure, ribozyme function and biological catalysis 
in general. 5 " 7 

Recent years have seen substantial progress in understanding of 
the structures and catalytic mechanisms of ribozyme systems. 8 " 14 
Among the best-studied catalytic RNAs are the self-cleaving 
ribozymes, including the hairpin and hammerhead motifs found 
in subviral RNAs from plant pathogenic viruses, the human 
hepatitis delta ribozyme, the VS ribozyme from Neurospora crassa 
and the ligand-activated glmS riboswitch. The hairpin ribozyme 
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from the satellite RNA of the tobacco ringspot virus is of par- 
ticular interest from a structural perspective because its active site 
is formed by the cation-driven interaction of two RNA internal 
loops to form a tightly integrated structure. In other words, an 
obligatory RNA tertiary structure transition, termed docking, 
precedes the chemical step of catalysis. In addition, full activity 
is observed when the physiologically relevant activating cation 
Mg 2t is replaced with the substitution-inert metal ion cobalt(III) 
hexamine [Co(NH 3 ) 6 3 *] 15 " 17 or with very high concentrations of 
monovalent salts in the complete absence of polyvalent cations, 18 
implying that metal ion cofactors do not play a direct role in 
reaction chemistry. The simplest interpretation of these results is 
that multivalent metal ions are required only for the formation 
of tertiary docked structure, but a role for metal ion cofactors in 
electrostatic stabilization, required conformational changes fol- 
lowing docking or other reaction steps has not been ruled out. 

The hairpin ribozyme exists naturally with loops A and B on 
adjoining arms of an RNA four-way junction, the presence of 
which greatly favors the formation of the docked state (Fig. l). l; " 
21 The ribozyme is also highly active in a two-way junction or 
"minimal" form, on which much of the biochemical work has 
been performed. 22,23 Finally, cleavage activity has been observed 
with internal loops A and B on separate molecules 24,25 (Fig. 1). In 
this junctionless form, which we refer to as the "trans-docking" 
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Figure 1. Hairpin ribozyme loop A and loop B constructs used in the 
present work shown with their context in the native, four-way junction 
ribozyme. Arrow indicates the site of self-cleavage. 



version of the ribozyme, the docking transition represents the for- 
mation of an RNA tertiary-structure interaction in isolation from 
hybridization or other effects. 

The hairpin ribozyme has been the subject of detailed 
analysis by a wide-range of biochemical and biophysical tech- 
niques, 10 ' 11,22,26 " 36 notably including extensive single-molecule 
analysis of the docking transition. 37 " 40 NMR solution structures 
of the isolated loop A and loop B constructs have been solved, 41 " 43 
and numerous X-ray crystallographic studies of the docked form 
have been performed. 11,44 " 52 Dramatic differences between the 
structures of the internal loops and of the complete ribozyme are 
observed, including the disruption of five of the seven nonca- 
nonical base pairs present in the solution structure of loop B, 44 ' 53 
implying extensive conformational rearrangements accompa- 
nying loop-loop interaction (Fig. 2). An intricate interdomain 
interface is formed that includes the formation of a ribose zipper, 
the creation of a pocket in which nucleotides from both loops 
envelop the base of U42 and the insertion of the loop A nucleo- 
tide G+l (immediately 3' to the cleavage site) into the stacked 
structure of loop B to form a base pair with C25. The reactive 
groups take up the catalytically relevant in-line attack geometry 
in the docked structure, consistent with docking being obliga- 
tory to the formation of the active state of the ribozyme. The 
mechanism of domain docking is thus clearly much more com- 
plicated than a rigid lock-and-key interaction and likely involves 
induced fit and/or conformational-capture processes on one or 
both loops. Low 3> values for the formation of individual interac- 
tions measured using single-molecule techniques imply an early 
transition state for docking, 54 but the molecular details of the 
process and its driving force are not well-understood. 

Biophysical studies of the hairpin ribozyme in the trans- 
docking form have the potential to allow new insight into the 
docking transition using direct measurements of bimolecular 
equilibrium and rate constants for binding and differential iso- 
tope labeling for spectroscopic purposes. Crystallography has 
confirmed that the interface in this format is essentially identi- 
cal to that in the native four-way junction. 48 Although a binding 
constant for the loop-loop interaction under conditions of high 
metal ion concentration has been reported (see below), detailed 



biophysical characterization of this form of the ribozyme thus 
far is minimal. 55 In particular, there have been no reports of the 
metal dependence of either docking or catalysis, or of the kinetic 
aspects of docking, for the trans-docking reaction. 

In this paper, we report equilibrium and kinetic analysis of 
the Co(NH 3 ) s 3t -driven intermolecular docking transition in 
the hairpin ribozyme. To identify ionic conditions under which 
stable tertiary structure forms, we use a novel circular dichroism 
(CD) difference assay to isolate the RNA-cation interactions that 
are responsible for driving intermolecular docking. Using surface 
plasmon resonance (SPR), we observe a relatively strong (sub- 
micromolar binding interaction between the two domains, 
accompanied by remarkably slow on- and off-rates that are con- 
sistent with the dramatic structural rearrangements accompa- 
nying docking. Finally, we use assays of self-cleavage activity to 
compare the ionic requirements for docking and catalysis and 
find that, despite being unable to participate directly in reaction 
chemistry or Lewis acid catalysis, Co(NH 3 ) 6 3+ plays additional 
roles in the catalytic mechanism of the hairpin beyond simply 
driving the docking transition. 

Results 

Intermolecular domain docking is driven by sub-millimolar 
concentrations of exchange-inert trivalent metal complex. CD 

spectra of nucleic acids are sensitive to overall structural fea- 
tures, including base stacking, and have found use in studies 
of RNA and DNA structure and folding. 56 " 58 The CD spectra 
of loops A and B of the hairpin ribozyme would be expected to 
be altered upon metal ion-driven domain docking, potentially 
providing a useful signature for the metal dependence of dock- 
ing. In a simple titration of a mixture of the two loops with 
Co(NH 3 ) s 3t , however, this effect would be obscured by signals 
arising from the extensive non-specific interactions of trivalent 
metal ions with the polyanionic RNA backbone. Indeed, CD 
spectra of isolated hairpin loops showed features consistent 
with A-form structure that increased in intensity upon equili- 
bration with sub-millimolar amounts of Co(NH 3 ) 6 3t (Fig. SI). 
We therefore devised a novel difference CD assay to isolate the 
cation-binding isotherm of only those metal ion(s) whose inter- 
action with the RNA is linked to domain docking. In this pro- 
cedure (Fig. 3), separate CD spectra at identical Co(NH 3 ) 6 3+ 
concentrations are acquired for three samples: Loop A alone, 
loop B alone and a 1:1 mixture of the two molecules. Samples 
must be carefully dialyzed against Co-containing buffer to 
ensure identical concentrations of free cation (see Materials 
and Methods). If no interaction between the two loops occurs, 
the CD signal for the 1:1 mixture [Ae (A+B) ] will be identical 
to the numerical sum of the signals for the two isolated loops 
(Ae A + Ae B ) (Fig. 3, upper panel). To the extent that the loops 
interact, however, the CD spectrum of the mixture will differ 
from the sum of the individual loops (Fig. 3, lower panel). The 
magnitude of the difference AAs, taken at the wavelength peak 
of nm, thus provides a direct readout of intermolecular domain 
docking. A plot of this parameter vs. [Co(NH 3 ) 6 3+ ] reveals an 
apparent cation-binding isotherm for only those cations whose 
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Figure 2. (A) Schematic of conformational changes upon formation of docked structure in the hairpin ribozyme, derived from NMR structures of the 
free loops 4143 and a 2.05 A crystal structure of the docked form in a junctionless construct. 48 Open ribose symbols are C3'-endo, shaded are C2'-endo 
and partly shaded (NMR structures only) display J-coupling constants characteristic of an equilibrium between the two conformations. Shading on 
the base moiety represents a syn conformation. The three major components of the interdomain interface are the G+1 :C25 Watson-Crick interdo- 
main base pair (red), a ribose zipper (green) and a complex pocket surrounding the extruded base U42 (blue). Base pair symbols are according to the 
convention of Leontis and Westhof, 83 in which filled symbols represent c/s base pairs, open symbols represent trans (including sheared purine:purine) 
and circles, squares and triangles represent Watson-Crick, Hoogsteen and sugar faces, respectively. Base pairs for the crystal structure were assigned 
using the RNAview annotation software. 84 Asterisks (orange) indicate observed Co(NH 3 ) 6 3+ ions, located in the major groove in all cases. The crystal- 
lography ion in loop B of the docked state is consistent with the location of Ca 2+ ions in the structure of the four-way junction ribozyme 44 and with 
lanthanide ion cleavage data. 85 The self-cleavage site (loop A), NOEs inconsistent with the major conformer (loop B) and in-line attack conformation 
(docked structure) are marked. (B) X-ray crystallographic structure of the junctionless hairpin ribozyme (PDB 20UE), 48 with interface elements colored 
to correspond with (A). Figure constructed with PyMOL. 86 



binding is thermodynamically linked to domain docking. The 
difference spectrum, unlike the CD signals of the individual 
samples, is unaffected by RNA-cation interactions that occur 
independently of domain docking. 

Representative data obtained under two ionic conditions are 
shown in Figure 4. In the absence of Co(NH 3 ) 6 3+ , domain dock- 
ing was not expected to occur, and the spectrum of the mixed 
sample was indeed essentially identical to the sum of the spectra 
of the isolated loops throughout the region of the ellipticity max- 
imum. Thus, AAe measurements at a particular [Co(NH 3 )/ + ] 
may be presumed to arise from the effects of metal ions only. In 
the presence of 250 |xM-free Co(NH 3 ) 6 3t , by contrast, a substan- 
tial shift between the two spectra was observed, leading to an 
easily quantifiable difference spectrum. 



Measurement of AAe 270 for our hairpin ribozyme constructs 
as a function of [Co(NH 3 ) 6 3+ ], although perturbed by significant 
data scatter, resulted in an unambiguous cation-binding isotherm 
with a [Co] 1/2 dock of 48.7 ± 35.0 (jlM (Fig. 5). No evidence for 
cooperative binding of multiple Co(NH 3 ) 6 3t ions was seen, and 
a Hill plot of the cation-binding data (Fig. 5, inset) had a slope 
of 0.87 ± 0.63. Comparable data for the G+1 A mutant of loop A 
showed no consistent ion-dependent changes in AAe 270 (Fig. SI), 
consistent with the incompetence of this mutant for docking and 
strongly favoring the interpretation that the effect observed for 
wild-type sequences arises from authentic domain docking. 

Interestingly, the Co(NH ) s 3+ requirements for trans-docking 
reported here are roughly comparable to those found by Burke 
and coworkers in a two-way-junction construct, which displayed 
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Figure 3. Schematic of the difference CD assay for the binding of ions linked to domain docking. In 
the absence of any docking, the observed CD spectrum of a mixture of loops A and B [Ae y solid line] 
will be the arithmetic sum of those for the individual loops (Ae A and Ae b , dashed lines) (top), whereas 
if docking occurs, a change in the ellipticity from that simple sum will be observed (shaded area, AAs) 
(bottom). Circles indicate multivalent cations bound to RNA, shaded in the case of ions whose binding 
is linked to domain docking. 



captured on the second channel as 
a reference. Figure 6A shows repre- 
sentative raw wild-type and mutant 
injection data. Both species were 
affected by a rapidly appearing bulk 
shift in refractive index, which we 
speculate arose from increased local 
concentration of high-mass cobalt 
complexes. A much slower, authentic 
macromolecular binding event was 
clearly present in channels contain- 
ing wild-type loop A but absent in 
the case of the non-docking G+1A 
mutant. Given the dramatic dif- 
ference in results upon this highly 
conservative sequence change, we 
interpret the slower event as repre- 
senting the specific formation of hair- 
pin ribozyme tertiary structure via 
intermolecular domain binding. 

A complete SPR analysis of 
loop A-loop B binding at 250 jjlM 



Co(NH, 



at which docking 



a [Co] yj dock of 17 ± 0.5 jJtM with n below l. 55 By contrast, the 
Lilley group reported data yielding a [Co], A dock of 0.24 jjlM for a 
four-way junction construct under somewhat different solution 
conditions (90 mM tris buffer pH 8.3). 59 In the latter case, a Hill 
coefficient of 1.75 was observed, indicating the presence of mul- 
tiple cooperatively binding ions, but this effect was attributed to a 
known ion-binding site in the four-way junction itself, a structure 
not present in our experiments. 60 

Slow association and dissociation kinetics in the trans- 
docking hairpin ribozyme. The above results demonstrate that 
intermolecular domain docking is driven by relatively modest 
concentrations of free Co(NH 3 ) 6 3+ . The kinetics and thermody- 
namics of trans-docking can thus be analyzed under saturating 
concentrations of trivalent metal ion. Fe-EDTA probing has pre- 
viously been used to investigate the trans-docking of loops A and 



B under high-[Co(NH 3 ) 6 3+ ] conditions, 55 resulting 



in a Kj°°v 



of 



4.8 ± 1.8 (jlM. No direct measurements of docking kinetics have 
been reported in the junctionless form of the ribozyme, although 
comparisons to off-rates observed in the single-molecule context 
are relevant (see below) . 

For analysis of the kinetics of hairpin ribozyme loop-loop 
interactions in the trans-docking format, we used a surface 
plasmon resonance (SPR) assay to monitor the formation of the 
bound complex in real time. Chemically synthesized 5 -biotinyl- 
ated wild-type loop A, rendered unreactive with a 2'-0-methyl 
modification at the nucleophilic group, was captured on one of 
the two SPR flow cells. Loop A with a G+1A mutation that ren- 
ders the molecule inactive for both docking and catalysis was 



reactions are driven essentially to 
completion (see above), is shown in 
Figure 6B. For this purpose, we plot- 
ted and fit the difference between SPR 
traces for wild-type loop A and G+1A 
mutant channels. This procedure yielded single-exponential data 
that could be interpreted directly in a l:l-binding model. We 
observed unusually slow on- and off-rates for macromolecular 
interactions, with mean association rate k a = (1.97 ± 0.29) x 10 3 
M" 1 s" 1 and dissociation rate k i = (7.1 ± 1.0) x 10" 4 s" 1 , resulting in 
an affinity of K d = 372 ± 101 nM. The slow kinetics of this bind- 
ing event, with an on-rate several orders of magnitude below the 
macromolecular diffusion limit, emphasizes the extensive con- 
formational changes required in each domain for the formation 
of the docked structure (Fig. 2). The observed k a is noticeably 
slower than previously studied RNA-RNA interactions includ- 
ing, for example, those for the kissing interactions between the 
HIV TAR RNA hairpin and either an RNA aptamer (9.2 x 10 4 
M" 1 s" 1 at 23°C) 61 or its complementary hairpin TAR* (2.1 x 10 5 
M" 1 s" 1 at 4°C). 62 Both of the latter processes require the unwind- 
ing of a duplex of at least five base pairs present in the unbound 
state for binding to occur. 

Our measured K A for the trans-docking construct is signif- 
icantly more favorable than that of 4.1 ± 1.8 |jlM reported by 
Burke and coworkers using Fe-EDTA footprinting to assay the 
formation of the ribozyme solvent-protected core. 55 Since core 
protection is assumed to correlate directly with the docking inter- 
action, this difference may arise from slight differences in ionic 
conditions or from minor variations in the construct, notably the 
use of a 2'-deoxy rather than 2'-0-methyl modification at the 
cleavage site in the footprinting work. Single-molecule kinetic 
fingerprinting results 63 have indicated a role for substituents at 
this group in modulating the loop interaction, and the precise 
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Figure 4. Comparison of CD spectra for an equimolar mixture of loops A and B [Ae (A+B) , solid line] with the arithmetic sum of spectra for the individual 
loops (Ae a + Ae b , dashed line). The difference between these traces (AAe, dotted line) is taken as a measure of the extent of ion-driven domain dock- 
ing. (A) Absence of Co(NH ) 3+ ; (B) 250 (jlM free Co(NH ) 3+ . 



effects of these modifications on on- and off-rates as well as 
on affinity are an interesting topic for further investigation. 

Increased metal requirements for self-cleavage com- 
pared with docking. The metal activation of hairpin ribo- 
zyme self-cleavage in two- and four-way junction constructs 
has been extensively studied (see Introduction). In situa- 
tions including the activation of a four-way junction con- 
struct with various metal ions, distinct ion dependencies for 
folding and cleavage have been reported. 20,59 The use of an 
intermolecular docking system in the present work allows 
a straightforward comparison of the metal dependence of 
domain docking, as determined above using CD, with that 
of ribozyme self-cleavage. These measurements take place 
in the absence of effects mediated by metal interactions 
with the junction. By considering the loop B construct as 
the enzyme and the loop A construct as the substrate, we 
determined the rate of ribozyme reaction k ohs under single- 
turnover conditions (excess of loop B over loop A) as a func- 
tion of Co(NH 3 ) 6 3+ concentration. In all cases, essentially 
complete cleavage of loop A was observed at long time points 
(Fig. 7), allowing us to consider the undocking of cleaved 
loop A from loop B and/or the dissociation of cleavage prod- 
ucts to be essentially irreversible. We thus neglected the 
effects of back-reaction in the kinetic analysis (see below) . 
We performed most of our measurements at 500 nM loop B, 
which is only slightly above the K A for the two loops determined 
above. However, when the 350 |jlM Co(NH 3 ) 6 3+ reaction was 
repeated at a loop B concentration of 3 |JlM, the resulting rate was 
identical within error to that observed at 500 nM loop B (data 
not shown), indicating that saturating concentrations of loop B 
had been reached. Consistent with this, a 10-fold lower apparent 
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Figure 5. Binding of ions linked to domain docking as measured by differ- 
ence CD analysis. Solid line, fit to equation 1. Inset, Hill plot of data through 
250 |jlM, showing a linear fit (solid line) and a linear fit with slope forced to 
unity (dashed line) for comparison. Data points represent the average of 
results for two independently prepared samples; error bars were derived 
via propagation from the reproducibility of measurements on individual 
samples. 



cation, k obs reached approximately 0.1 min" 1 (extrapolated value 
0.113 ± 0.027 min" 1 ), consistent with limiting values reached in 
previous studies of this format. 25,55 No cleavage was observed for 
the control G+1A mutant of loop A (data not shown). 

Our measurements showed a strongly sigmoidal dependence 
of k obs on [Co(NH 3 ) 6 3+ ], implying cooperative binding of mul- 
tiple ions necessary for cleavage (Fig. 8). Fitting with Equation 2 



K vs. docking K d was observed by Burke and co-workers, 55 pos- gave a [Co].. m of 168 ± 51 U.M. This value is significantly 



sibly attributable to the difference in chemical structure at the 
nucleophilic oxygen (see above). At saturating concentrations of 



higher than the [Co] 1/;! dock observed above for interdomain dock- 
ing as assayed by CD, implying that metal ions play a role in 
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Figure 6. (A) Unreferenced surface plasmon resonance data acquired at 250 |j,M Co(NH 3 ) 6 3+ and 450 nM loop B. Upper trace, wild-type loop A; lower 
trace, G+1A mutant loop A immobilized on the same chip. (B) Representative fully referenced SPR analysis of domain docking at 250 (jlM Co(NH 3 ) 6 3+ . 
Data shown are differences between a channel with immobilized wild-type loop A and one with immobilized G+1A mutant (compare left panel) and, 
thus, reflect only specific docking interactions. Solid lines are global fits to a 1:1 interaction model incorporating bulk refractive-index shift. Concentra- 
tions of loop B used: 150 nM, 300 nM, 450 nM, 600 nM, 750 nM. 
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hairpin catalysis in this format beyond simply driving the trans- 
docking interaction. Consistent with this interpretation, the cor- 
responding Hill plot exhibited linear behavior throughout the 
measured range with slope n of 1.87 ± 0.15, contrasting with the 
non-cooperative behavior observed for docking. 



Interestingly, the [Co] yj for catalysis measured here is 
roughly comparable to the value of 0.1 mM for a two-way junc- 
tion construct (in 50 mM tris buffer pH 7.5) reported by Fedor 
and co-workers. 16 This result may imply that, although the 
ion requirements for domain docking differ greatly according 
to the presence and nature of the interdomain junction, 59 the 
binding of the final ion(s) necessary to effect catalysis proceeds 
with relatively little influence of the junction once the docked 
region has been formed. The role of the junction sequences in 
the hairpin ribozyme would thus contrast somewhat with that 
of the tertiary-stabilizing elements in the full-length hammer- 
head ribozyme or the kissing-loop tertiary interactions in the VS 
ribozyme, both of which act to dramatically decrease the ionic 
requirements for catalysis as well as folding. 64,65 Substantial 
further work would be necessary to explore this possibility 
fully, especially given the slightly different solution conditions 
between the data sets in different constructs. 

Discussion 

We report a detailed analysis of the thermodynamics and kinet- 
ics of intermolecular domain docking in the hairpin ribozyme 
activated by the exchange-inert metal complex Co(NH 3 ) s 3+ . We 
found high-affinity domain docking to be driven by sub-milli- 
molar concentrations of free trivalent cation, with no detectable 
binding for a mutant previously determined to be deficient in 
docking, emphasizing the strong similarities between the docked 
states formed in native, minimal and junctionless ribozymes. 46,48 
A comparison between the activation curves for stable docking and 
RNA self-cleavage (Table 1) revealed a greater ionic requirement 
for the latter, indicating a role in catalysis for ions binding more 
weakly than those directly coupled to tertiary structure formation. 
Similarly, the self-cleaving hammerhead ribozyme, in which metal 



Figure 7. Cleavage kinetics of the trans-docking hairpin ribozyme at 250 
jjiM Co(NH 3 ) 6 3+ . Top, 20% polyacrylamide gel analysis of the cleavage of 
3'- 32 P-labeled loop A RNA, demonstrating the lack of starting material 
present at long time points. S, uncleaved loop A substrate (26 nucleo- 
tides); P, 3' cleavage product (21 nucleotides). Bottom, densitometric 
scan of the gel data fit to a single rising exponential. For this data set, the 
extrapolated cleavage at infinite time is 86.1%. 
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Figure 8. (A) Metal-dependent kinetics of the hairpin cleavage reaction in trans-docking format. Solid line is the fit to Equation 2, whereas the dashed 
line is the best-fit CD analysis of docking from Figure 5 scaled to the same limiting y-value at [Co]->=° as the kinetic data. (B) Kinetic data shown as a 
Hill plot. (C) Kinetic scheme for the reaction in trans-docking form indicating the observed irreversibility of undocking and/or dissociation of cleavage 
products (see text). Rate constants k a and k 6 correspond to those measured using SPR in this work. 



ions appear to play a more direct role in catalysis than in the hair- 
pin, has shown a distinction between ionic requirements for ter- 
tiary structure formation and self-cleavage in several contexts, 66,67 
and the identity of the divalent metal used was found to strongly 
affect the efficiency of catalysis but not tertiary structure formation 
in the hammerhead. 67 In the four-way junction form of the hairpin 
ribozyme, multiple phases of ion-driven structure formation and 
distinct ionic requirements for docking and catalysis have been 
previously reported for Mg 2t -activated cleavage. 20,68 The present 
results refine this picture by ruling out models that explain the dif- 
ferential metal dependence by invoking ion interactions with junc- 
tion sequences, conformational changes driven by inner-sphere ion 
interactions with RNA functional groups or metal ion participa- 
tion in catalysis via inner-sphere effects. It is also interesting to 
note that Mn(II), although it does not support hairpin ribozyme 
self-cleavage, nevertheless leads to the formation of global tertiary 
structure indistinguishable from that supported by Mg 2t in the 
four-way-junction form of the ribozyme, 5 ' consistent with domain 
interaction being a necessary but not sufficient condition for cat- 
ion support of catalysis in this system. The catalytic role played by 
ions not linked to domain docking is not known, but may involve 
an ion-induced conformational change of the docked structure or 
electrostatic participation in catalysis by ions not directly coordi- 
nated to the reactive groups. 6,59 

Although the internal equilibrium between cleavage and 
ligation of docked hairpin ribozymes is of order of magnitude 



unity under many conditions, and indeed often favors liga- 
tion, 65 we observed essentially complete cleavage of loop A at 
long time points in all of our kinetic assays (Fig. 7). We inter- 
pret this observation as arising from rapid, irreversible undocking 
and/or dissociation of cleavage products. Consistent with this, 
undocking in the two-way junction form was found via kinetic 
fingerprinting to be accelerated by two orders of magnitude in 
the cleaved (0.043 sec' 1 ) vs. the uncleaved species. 63 Likewise, 
undocking in the four-way-junction form is dramatically accel- 
erated upon substrate cleavage to a rate much greater than that 
of religation. 70 Unless the cleavage product is much more stable 
to undocking in the absence than in the presence of the highly- 
structured junction, which we consider unlikely on structural 
grounds, the rate of product release in our experiments can thus 
be presumed to be much greater than that of religation across the 
full-range of [Co(NH 3 ) 6 3+ ] values used, allowing back-reaction 
to be neglected. Thus, although the present results do not iden- 
tify the elementary reaction step responsible for the observed ion 
dependence of £ obs , the comparisons between k ohs and extent of 
docking as measured by CD may be interpreted in straightfor- 
ward fashion. The fundamental observation that at intermedi- 
ate [Co(NH 3 ) 6 3t ] (e.g., 75 (JtM to 100 |xM in Fig. 8) the extent 
of docking has reached a much higher fraction of its saturating 
value (8 = 1) than has the cleavage rate strongly supports our 
qualitative conclusion that cations play a role in the support of 
catalysis beyond that of simply driving domain docking. 
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Table 1. Summary comparison of metal dependence of domain docking 
and catalysis 





[Co]% (|j,M) 


n 


Docking (difference CD assay) 


48.7 ± 35.0 


0.87 ± 0.63 


Catalysis (self-cleavage assay) 


167.9 ±51.1 


1.87 ±0.15 



It should be noted that phenomenological cation-binding 
isotherms such as those reported here do not uniquely demon- 
strate that individual ion(s) bind to well-defined site(s) on the 
nucleic acid structure. Especially under low-salt conditions, elec- 
trostatic accumulation of cations can occur in non-mass-action 
driven fashion, in which case [Co]^ and n are phenomenological 
descriptors rather than true molecular representations of K i and 
cooperativity, respectively. 2,3,71 In the isolated loop A, ion titra- 
tion at concentrations leading to a substantial change in AAe at 
270 nm (Fig. SI) did not significantly alter the aromatic-ribose 
fingerprint region of a 1 H- 13 C HSQC NMR spectrum (Fig. S2). 
This observation is consistent with the lack of change in NMR 
data for a somewhat different loop A construct upon addition of 
3 mM Mg 2 *, 41 as well as with the lack of paramagnetic broaden- 
ing observed upon titration of loop A with Mn(II), 43 although 
both of the latter results were obtained at higher concentrations 
of monovalent salt than were used in this work. Regardless of 
the physical interpretation of the binding isotherm, however, 
the comparison of ionic requirements for docking and catalysis 
(Fig. 8, Table 1) remains valid. 

Using direct on-rate measurements via SPR, we found a 
remarkably slow association rate between the two domains con- 
sidering the relatively tight binding constant. Based on the sub- 
micromolar binding constant, we assume this association rate 
to correspond directly to the formation of the docked structure 
seen crystallographically and, thus, to probe the same process as 
the CD experiments. Our construct lacks the four-way junction 
that has been shown to be the source of the approximated but 
not docked state visible in some single-molecule FRET experi- 
ments. 3 ' A value of k five to six orders of magnitude below the 

on o 

diffusion limit implies, in the simplest interpretation, that only 
a very small fraction of encounters between loops A and B are 
ultimately productive for the formation of the loop-loop interac- 
tion. 72 This observation may be understood in light of the fact 
that both binding partners have substantial albeit dynamic struc- 
ture in the undocked state (Fig. 2). In this situation, high-affinity 
binding is maintained via an extensive intermolecular interface 
leading to a very slow off-rate. 

Although these results represent, to our knowledge, the first 
direct measurement of second-order rate constants for domain 
association in the hairpin ribozyme, it is interesting to compare the 
off-rate measured here with undocking rate constants (^ undock ) pre- 
viously observed in unimolecular (two-way junction) constructs. 
Single-molecule fluorescence resonance energy transfer (FRET) 
experiments have found four distinct docked states with k , , 

i undock 

varying from 0.005 sec"'-3 sec" 1 for 2'-0-methyl modified ribo- 
zymes in 12 mM Mg 2 *. 37,40 Subsequent kinetic fingerprinting stud- 
ies determined that k , , for the most stable population decreased 

undock r r 

to 0.00045 sec" 1 for a native, cleavable (2'-OH at A-l) substrate, 



emphasizing the effect of modifications at the cleavage site on the 
nature of the docked state. 63 The 5 -fold lower off-rate measured 
here compared with £ UIldock for a two-way junction construct with 
2'-0-methyl at A-l may most likely be attributed to the greater 
charge on Co(NH 3 ) 6 3+ compared with Mg 2+ . With this important 
caveat, we see no evidence for specific stabilization of the docked 
state by the two-way junction, consistent with the view of this 
structure as indistinguishable from a flexible linker. 73,74 

The measured K d of 372 ± 101 nM corresponds to a AG° of 
-8.91 ± 0.19 kcal/mol at 30°C. This value is somewhat more nega- 
tive than that typically observed in the single-molecule experiments 
described above, again perhaps due to different ionic conditions. 
Thermodynamic measurements on RNA tertiary structure for- 
mation are relatively rare, but have recently been reported for the 
GNRA tetraloop-tetraloop receptor complex using isothermal 
titration calorimetry (ITC) 75 and single-molecule 76,77 techniques. 
Calorimetrically determined AG° values in this system vary from 
-10.0 kcal/mol— 9.6 kcal/mol between 40°C and 47.5°C for a biva- 
lent construct, suggesting a AG° of -5 kcal/mol to -6 kcal/mol for 
a single tetraloop-tetraloop receptor interaction at 30°C. 75 Also, for 
the GNRA tetraloop-tetraloop receptor complex, single-molecule 
techniques measured a k off of more than 3 sec" 1 at saturating Mg 2t , 
several orders of magnitude higher than that reported here for the 
hairpin, 78,7 ' emphasizing the tightly interconnected structure of 
the docked form of the ribozyme. 

The very slow on-rate in our experiments could be explained 
by a model in which both isolated loops dynamically sample con- 
formations resembling their docked states and only collisions in 
which each loop is simultaneously in such a state are productive 
for binding. Such a "double conformational capture" mechanism 
is not applicable in the arguably more common case in which 
only a single binding partner has pre-existing structure that must 
be disrupted for binding to occur, consistent with the relative 
rarity of slow binding kinetics such as observed in the hairpin. 
The observation of individual NOEs in isolated loop B that are 
consistent with the docked structure supports this possibility. 42,44 
Taking advantage of the biophysical information reported here 
and the spectroscopic flexibility afforded by the trans-docking 
format, we are currently engaged in a systematic study of the 
dynamic properties of the isolated loops in order to explore the 
double conformational capture hypothesis. 

Materials and Methods 

RNA preparation. The following RNA sequences 
were used: Loop A (wild-type), 5-GCG 
CAmGUCCUCGUAAGAGAGAAGCGC-3'; Loop A (G+1A), 
5'- GCG CAmAUCCUCGUAAGAGAGAAGCGC-3'; Loop 
B, 5'-GCGAGAGAAACACACGACGAAAGUCGUGGU 
ACAUUACCUCGC-3'. Am indicates a 2'-0-methyl-adenosine 
modification used to preclude cleavage at the scissile phos- 
phate. This modification was omitted in loop A sequences used 
for ribozyme cleavage assays. Wild-type and mutant loop A 
sequences (26 nt) were chemically synthesized by the W.M. Keck 
Foundation Biotechnology Resource laboratory at Yale University, 
deprotected using tetrabutylammonium fluoride (TBAF) solution 
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in tetrahydrofuran (THF) for 6 h at room temperature and 
desalted on Sephadex G-25 (Sigma-Aldrich). Loop B (42 nt) was 
synthesized by in vitro transcription from a synthetic DNA tem- 
plate with T7 RNA polymerase 80 and desalted via ethanol precipi- 
tation. Desalted RNAs were purified on a Superdex 75 (26/60) 
column (GE Healthcare) in which the eluent was 10 mM phos- 
phate buffer (pH 6.5) and 100 mM NaCl at a flow rate of 3.0 mU 
min. 81,82 RNA-containing fractions were combined, concentrated 
and desalted against RNase-free, double-deionized water using 
Amicon Ultra 15 centrifugal filtration units, dialyzed extensively 
against 20 mM HEPES, 20 jjlM Na 2 EDTA at pH 7.5 (HPLS 
buffer) and finally dialyzed against HPLS containing the desired 
concentration of Co(NH 3 ) 6 3+ . Dialysis was performed on a stirred 
Spectra/Por microdialyzer with a regenerated cellulose membrane 
of molecular-weight cutoff 2000. RNAs were renatured by heating 
to 65°C for 5 min followed by quickly cooling on ice for 10 min. 
For immobilization on optical chips, 5 -bio tin-tagged RNAs 
were synthesized by Dharmacon, Inc. and deprotected accord- 
ing to the manufacturer's protocol. RNAs were then lyophilized 
and exchanged into RNase-free, double-deionized water using 
Amicon Ultra 4 centrifugal filtration units, analyzed for purity on 
a denaturing polyacrylamide gel, lyophilized to dryness and stored 
at -20°C. RNA was quantitated using extinction coefficients at 
260 nm determined by comparison of absorbance before and after 
PI nuclease digestion using the known absorbance values of the 
mononucleotides. 56 Final values were 214,000 M" 1 cm" 1 for wild- 
type loop A, 232,000 M" 1 cm' 1 for G+1A and 354,000 M' 1 cm 1 for 
loop B. For CD samples containing 1:1 mixtures of two compo- 
nents, a final dialysis of at least 1 h against the desired concentra- 
tion of Co(NH ) 6 3+ was performed following mixing of the two 
components and the final concentration was determined using the 
sum of the individual extinction coefficients (i.e., changes in e 
upon formation of docked structure were neglected) . 

Circular dichroism spectroscopy. Circular dichroism spec- 
tra were obtained with a Chirascan spectropolarimeter (Applied 
Photophysics). Experiments were performed in duplicate using 
independently dialyzed samples in HPLS buffer at 30°C and 
2-3 ixM RNA. Data were collected from 200-400 nm, aver- 
aged over four scans and converted to molar ellipticity (Ae) using 
extinction coefficients given above. Control scans of Co(NH 3 ) s 3+ - 
containing buffer samples showed negligible CD signal, justifying 
the neglect of contributions from free Co(NH ) 6 3t . The binding of 
docking-related ions was identified according to the molar elliptic- 
ity difference (AAe) at 270 nm between Ae for a sample of mixed 
loops A and B and the sum of Ae values for individual loop A and 
loop B samples. The mean deviation between duplicate determina- 
tions of AAe at a single metal concentration was 12.7%. The data 
were fit using nonlinear regression (Igor Pro 5.04B) to where [L] f is 

[L], 



AAs = AAe„ 



[Co]?;? +[L] f 



(1) 



the concentration of free Co(NH 3 ) 6 3+ as fixed by sample dialysis. 
Hill analysis was performed by defining 0 as the ratio between 
AAe at a given [Co(NH 3 ) 6 3t ] and the AAe max extracted from the 
binding-isotherm curve fitting. As the measured value for AAe at 
500 |xM [Co(III)] exceeded the extrapolated AAe max , that data 



point was perforce excluded from the Hill analysis. Errors are as 
reported by the fitting routine. 

Surface plasmon resonance. Binding kinetics studies were 
done using a Reichert Technologies SR7000DC two-channel 
Surface Plasmon Resonance System at 30°C. The sensor chip 
used was a mixed self-assembled monolayer [10% COOH-(EG) 6 - 
Alkanethiol, 90% HO-(EG) 3 -Alkanethiol] surface (Reichert) to 
which streptavidin (NEB) was covalently attached using standard 
EDC/NHS coupling. Biotinylated wild-type and non-dockable 
control (G+1A) loop A were immobilized on separate channels 
in HPLS buffer to typical densities between 100-125 ixRIU and 
subsequently equilibrated in HPLS with 250 ixM Co(NH 3 ) 6 3+ . 
Loop B samples were dialyzed for 1 h in HPLS buffer with 
250 jjlM Co(NH 3 ) 6 3t and diluted to the desired concentration 
in the same buffer for the kinetic experiments. SPR experiments 
were performed with a 15 min association phase and 10 min dis- 
sociation phase using a flow rate of 10 uX/min. Test scans indi- 
cated that results were independent of flow rate in this range. 
Regeneration of the baseline was achieved by injecting 1 M NaCl 
and ddH 2 0 for 10 min each. The experiment was repeated on a 
total of three independently prepared chip surfaces with two to 
three injections on each surface. 

Data were analyzed using Scrubber 2 (BioLogic Software) 
and fit globally for each series of injections using a simple 1:1 
Langmuir interaction model including bulk refractive-index 
shift. The signal from the reference channel [Loop A (G+1A)] 
was subtracted from the signal from the channel with WT Loop 
A to yield the specific docking. The average of multiple blank 
(buffer) injections was then subtracted from each curve. Rate and 
equilibrium constants are thus reported as mean and standard 
deviation of global fits to eight independent data sets. 

Ribozyme kinetics. Loop A (-20 pmol) was labeled at the 3' 
end with (5'- 32 P)pCp (Perkin-Elmer Life Sciences) and T4 RNA 
ligase (New England Biolabs Inc.) using standard methods and 
purified on 20% polyacrylamide gels. Cleavage reactions were 
performed in HPLS buffer at 30°C and were initiated by addi- 
tion of Co(NH 3 ) 6 3t to the desired total concentration. A typical 
80 |xL cleavage reaction contained ca. 5,000 cpm of labeled loop 
A and 500 nM loop B. For each time point, a 10 |xL aliquot was 
removed, quenched in an equal volume of gel-loading buffer (87% 
formaldehyde, 25 mM Na 2 EDTA, 0.02% each xylene cyanol and 
bromophenol blue) at 95°C for 2 min and stored on ice. The cleav- 
age reactions were analyzed on 20% polyacrylamide gels and visu- 
alized on a Molecular Dynamics phosphorimager (GE Healthcare) 
with ImageQuant software (Molecular Dynamics). All experi- 
ments were repeated two or more times. The rate of hairpin ribo- 
zyme cleavage at each concentration of Co(NH 3 ) s 3+ was calculated 
by fitting percent cleavage product present to a three-parameter 
rising single exponential equation (Igor Pro 5.04B). The extrapo- 
lated extent of cleavage at completion of reaction varied between 
75-98%, with all but one concentration (90 (jlM) exceeding 85%. 
The resulting plot of k obs vs. [Co(NH 3 ) 6 3t ] exhibited clear sigmoi- 
dicity (see Fig. 8) and was therefore fit to a cooperative binding 
equation: 

[L]" 



*obs °^( [Co ]^)» +[L] » 



(2) 
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where [L] is the concentration of Co(NH 3 ) 6 3t and n is a Hill 
coefficient (Igor Pro 5.04B). The reported Hill coefficient of 1.87 
± 0.15 was extracted from the slope of a linear Hill plot of the 
data and was identical within error to that obtained from direct 
fitting to Equation 2 (1.86 ± 0.55). Errors in [Co] v4 chem and n are 
as reported by the fitting routines. 
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